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Abstract 

We answer a question of Laszlo Babai concerning the abelian sandpile model. Given 
a graph, the model yields a finite abelian group of recurrent configurations which is 
closely related to the combinatorial Laplacian of the graph. We explicitly describe the 
group elements and operations in the case of thick trees with loops — that is, graphs 
which are obtained from trees by setting arbitrary edge multiplicities and adding loops 
at vertices. We do this both concretely (by describing the so-called recurrent and 
identity configurations) and abstractly (by computing the group's abstract structure), 
and define maps identifying the two. 

1 Introduction 

The abelian sandpile model (ASM) was defined in [8] for any graph, generalizing the case 
of a grid from [3]. This model is a prime example of self-organized criticality [3] which has 
transformed the understanding of how complexity arises in nature (see [18] for a summary 
of applications to modelling earthquakes, forest fires, landslides, etc.) The ASM is related 
to the chip- firing game, introduced by Spencer [IB] , and later modified by Bjorner, Lovasz, 
and Shor [5]. For general references on the abelian sandpile model, see e.g. [HE]. 

We will only consider undirected graphs. An (undirected) abelian sandpile model begins 
with any connected undirected graph (with arbitrary edge multiplicities) with vertex set 
V and edge set E, together with a distinguished vertex s e 7 called the sink. Then, a 
configuration of the abelian sandpile model is an element of Z>° , Vq : = which assigns 

a nonnegative integer to each vertex except the sink, considered the number of grains of sand 
at that vertex. A configuration is stable if the number of grains of sand on each vertex is 
less than the degree of the vertex (otherwise it is unstable). Then, one obtains a dynamical 
system on the space of configurations, where at each time step, each unstable vertex sends 
one grain of sand along each incident edge. Sending grains of sand to the sink merely reduces 
the total number of grains. Because every vertex has a path to the sink, it turns out that 
this yields a well-defined stabilization of any configuration in finite time. The system can 
remain in motion by continually adding grains (e.g., at random vertices). 
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One basic question about any dynamical system is to find "recurrent configurations." 
For sandpiles, Dhar [8\ defined these to be stable configurations that can be obtained from 
any configuration by adding grains and stabilizing. Dhar showed that they form an abelian 
group (under adding and stabilizing), called "the sandpile group," whose abstract group 
structure is Z n /AZ n , where A is the matrix obtained from the Laplacian of the graph by 
deleting the row and column corresponding to the sink. 

However, the concrete structure (in terms of configuration on the graph) of these recurrent 
configurations turns out to be quite complicated and interesting. For example, Creutz's paper 
[6] displays intricate fractal patterns of the identity element of this group in the case of a 
square grid, which was further studied in [?l [13]. 

Laszlo Babai suggested that the authors study the sandpile group for "thick" graphs: 
graphs which may have large edge multiplicities, but whose combinatorial structure is simple. 
In this note, we consider thick trees with loops: graphs which become trees when all positive 
edge multiplicities are dropped to one, and all loops are discarded. We concretely describe 
the group elements and operations in these cases, compute the abstract group structure, and 
relate the two using recursive formulas. To our knowledge, this is the first result giving the 
complete, explicit structure of the sandpile group in any class of undirected graphs. In [TU] 
(the grid) and [21CI21EL7] (other graphs), information about the abstract group structure for 
more complicated graphs was given. 

The note is organized as follows: in Section [2] we recall the necessary background. In 
Section [3] we state our main theorem and its corollaries. In Section [4] we prove the theorem. 

2 The abelian sandpile model 

We review the theory of the (undirected) abelian sandpile model (ASM), following pQ. 

Let A be a finite connected undirected graph with arbitrary edge multiplicities, allowing 
loops. Precisely, A = (V,E,g), where V is the set of vertices, E is the set of edges, and 
g : E — > is the incidence map. Here is the set of one- and two-element subsets of 
V (i.e., unordered pairs of vertices). We have g(e) = {i,j} where i, j are the endpoints of e. 
We require that A be connected, i.e., there is a path of edges connecting any two vertices. 

Definition 2.1. An (undirected) ambient space is a collection (A, s) where A is an undi- 
rected, connected graph with vertex set V and edge set E (allowing multiple edges and 
loops), and s G V is a distinguished vertex, called the sink. All other vertices are called 
ordinary. Let Vq :—V\ {s} be the set of ordinary vertices. 

Notation 2.2. We will use bold letters (e.g. u) to denote vectors, usually in Z>° . 

Definition 2.3. A configuration of the ASM associated to (A, s) is a vector u G Z>° , i.e., 
an assignment of a nonnegative integer Ui to each ordinary vertex % (called its height). 

Definition 2.4. For any i,j G V, define the edge multiplicity to be e^j := #{a G E : 
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We may now define the notion of stability of a configuration: 

Definition 2.5. Given a configuration u, a vertex i e Vq is called stable if Ui < deg(z). 
Otherwise i is called unstable. A configuration is called stable if all vertices are stable. 

Definition 2.6. If % £ Vo is unstable for u, then define the state Ti(u)i by 7"j(u)j := ttj — 
deg(i) + e^j and Ti(u)j := Uj + e^j for j ^ i. Passing from u to Tj(u) is called "toppling the 
vertex i." 

Given any unstable configuration, one may continually perform such topplings n until a 
stable configuration results, which we call the stabilization. 

Theorem 2.7. |2l (cf. also JZV) Given any ambient space (A,s) and any configuration 
u of its ASM, then there exists a unique stable configuration o~(u), which satisfies o~{u) = 
T h T i2 ■ ■ ■ T i m ( u ) f or some sequence of topplings r h , r i2 , . . . , T im . 

(The theorem rests on an application of the Diamond Lemma, since if two vertices are 
unstable at once, one may topple them in either order with the same result.) 

Definition 2.8. For each configuration u, we denote by o~(u) its unique stabilization. 

With this in mind, one may conclude that the set of stable configurations form a monoid 
under vertexwise addition: 

Definition 2.9. [8] Let M be the set of stable configurations. Define the operation © on 
M by u (B v := a(u + v), where + is vector addition. 

Lemma 2.10. JZ|/ The operation © is commutative and associative, making (M, ©) a com- 
mutative monoid, called the sandpile monoid. 

Even the abstract structure of the sandpile monoid (without discussing concrete configu- 
rations) can be quite complicated. However, Dhar [8] noticed that the subset of "recurrent" 
configurations has a much simpler abstract structure: 

Definition 2.11. A configuration u is called recurrent if, for all configurations v, there 
exists a configuration w such that v © w = u. 

Proposition 2.12. J2J/ The set of recurrent configurations, G, is an abelian group under © ; 
called the sandpile group. 

(Actually, as noticed in p], the above is true for any finite commutative monoid, if one 
replaces G by the unique minimal ideal.) In [8] (elaborated in [TP] ) there is a surprisingly 
simple formula for the abstract structure of this group: 

Definition 2.13. The toppling matrix A = (Ajj)j je y is given by 




(2.0.1) 
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Note that any two configurations which have the same stabilization are related by an 
element in the row span of A; also, adding a large enough multiple of the sum of the rows 
Aj must give a unique stabilization. Using this, one has the 

Definition 2.14. Let A := (Aj) ig y be the lattice spanned by the rows Aj of the toppling 
matrix. 

Theorem 2.15. |2l U\/ The natural map Z>° -» G, c i— > a(c + e A ) descends to an isomor- 
phism Z y °/A = G. In particular, = det(A). 

(Note: G = Z Vo /A and = det(A) are due to [8j; the map c ^ a(c + e A ) was 

pointed out in [I].) 

Dhar found an algorithm to test for recurrence, of linear time in the combinatorial size 
of A (the number of edges counted without multiplicity plus the number of vertices): 

Definition 2.16. Let (3 G Z^° be given by (3 = J2 ie v A i- That is ' A = e M- 
We consider adding (3 to be "toppling of the sink" . One has the 

Algorithm 2.17 (Burning algorithm). [HI [15] To test whether a stable configuration u is 
recurrent, stabilize u + (3. Each vertex can topple at most once in the stabilization process, 
and u is recurrent iff every vertex topples (equivalently, iff u © (3 = u). 

The term "burning algorithm" (with which each vertex that topples in the process is 
considered "burned") was introduced in [14] . The algorithm above was discovered by [8], 
and the proof follows from [8] and [T5] . 

3 The sandpile group for thick trees 

Definition 3.1. Following L. Babai, we call an undirected graph T a thick tree if the 
underlying simple graph T' obtained by reducing all nonzero edge multiplicities to 1 is a tree 
(i.e., a connected graph without cycles). A vertex i is called a leaf if deg T / i = 1 and i ^ s. 

Definition 3.2. An undirected graph T is called a thick tree with loops if, after removing 
all loops (dropping e^j to zero for all i 6 V), the obtained graph is a thick tree. 

Recall that a graph is a tree iff, for any two vertices, there exists a unique path from one 
to the other. Let i, j G V. We define a partial ordering y on V, where i >z j if j is on the 
unique path of T' from i to the sink. By definition, this means that j >z s for all j G V. 
Also, i y j if the condition i ^ j is added. Vertex p(j) is the parent of j if p(j) y j and 
p(j) is adjacent to j. 

For a thick tree with loops, T, we define y using the underlying tree T' . 

We have the following main result, whose proof is the content of the next section. 

Notation 3.3. For any vertex i G Vq, let Si be the elementary vector given by {Si)j = 5ij = 1 
if i = j and otherwise. 
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Theorem 3.4. The sandpile group (G C Z>° , ©) for a thick tree T with loops can be char- 
acterized as follows (recall Notation \3.3\) : 



(i) G = {ue Z^ : Vj G V , deg(j) > Uj > deg(j) - e jAj) }. 

(li) The map <p : G -»• Y[ j&Vo Z/e iiP(i ), u h-> ^ 6 y 5 j (E;^ M * ( mod e i, P (i))) ^ « ^ ro ^ 
isomorphism. In particular, the abstract group structure does not depend on which 
vertex is chosen to be the sink. 

(Hi) The map cf)^ 1 may be expressed recursively as follows: For any j G Vq, 



Vj, if j is a leaf, 

fj [ Vj - ^ , otherwise, (3.0.2) 



where for any m G 1>/p, we let m G {0, 1, ... ,p — 1} be its preferred representative, 
and for any j G V , fj : Z -> {deg(j) - e iiP(i) , . . . , deg(j) - 1} owen (uniquely) by 
fj(m) = m (mod e ijP(i )). 

Corollary 3.5. Lei u fre any configuration. Then the associated recurrent configuration 
a{u + e T ) is given recursively by 

a (u + e T )- = [ Uj ~ \- u il e i*®\ e is*i) %s a lea f (3 3) 

[Uj + Xje jjP (j) + J2 iyj (ui - a(u + e T )i) otherwise, 

where Xj = [(deg(j) — Uj — J2i^j( u i ~ °~( u + eT )«) — l) / e i,p(j)J ■ If °~{ u ) ^ s already recurrent, 
we can delete the e T 's. 

Proof. Simply apply the map (j) and then the map we then expand fj. □ 
Corollary 3.6. The identity configuration e T is given by 

T= |0 if j is a leaf, (3.0.4) 

{ VjVpCj) ~ z2i H e i otherwise, 

where Xj = [(deg(j) + J2 iyj ef ~ l )hj, P {i)\- 

Proof. Specialize Corollary 13.51 to the case of u = 0. (Alternatively, specialize (13.0.21) to the 
case of v = and expand fj.) □ 



4 Proof of Theorem 1374 

We divide Theorem 13 . 41 into Propositions 14.11 14.31 and 14.41 corresponding to (refinements of) 
parts (i), (ii), and (iii), respectively. 
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Proposition 4.1. In a thick tree T with loops, the following statements are equivalent: 

(i) Configuration u is recurrent. 

(ii) For all j G V , deg(j) > uj > deg(j) - e jiP(j) . 

Proof. Take any stable configuration u. We apply Algorithm 12. 171 This means that we try 
to stabilize u + f3, and test whether or not every vertex topples. Since (3 adds sand only to 
the vertices adjacent to the sink, one easily sees that no vertex can topple before its parent. 
We can view adding (3 as having "the sink topple," and thus inductively, a configuration is 
recurrent iff every vertex can topple once its parent topples. We conclude that every vertex 
fires iff, for every vertex j, the amount of chips needed to topple (deg(j) — Uj) is less than 
or equal to the amount of chips it inductively receives from its parent (e^)). □ 

Lemma 4.2. The lattice A can be characterized as 

Proof. First, we show the inclusion C. It suffices to show that, if v = A^ for any k 6 V , 
then e jtP (j) \ Y,i>j v i for an Y 3 e v o- 



For k = p(j), one has 



J2 Vi = - e i,p(i) = ( mod e i>f(i))- ( 4 - - 6 ) 



For k = j, one has 

^2 v i = e i, P (i) = ( mod HpU))- ( 4 - - 7 ) 



For all k eV \ {j,p(j)}, one has 



5>i = 0. (4.0.8) 

The result is proved. 

Next, we show the inclusion D. For any j G Vq, consider the sum 

E A « = - e j,PU)( s p(j) - S j)- ( 4 - - 9 ) 

i>j 

A simple inductive argument on Vo shows that these elements span the desired space W = 
{v : Sj }P (j) | J2iyj v i^3 e ^o}- To do this, note that if j is a leaf, Ajj = ej )P ^, and one 
may express all elements of If as a sum of such Aj's and vectors which are zero in the leaf 
components. □ 

As a consequence, we deduce the 
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Proposition 4.3. One has an isomorphism : Z/A = Y[jev ^/ e j^pU) w hich descends from 
the map 

: Z§j -»• Y[ Z / e j,p(j)i ( t ) ( u )j = ^2 u i ( mod e j, P (j))- (4.0.10) 

Hence, the restriction <p\ G : G — > YljeV ^/ e j,pU) ^ s a 9 rou P isomorphism, and one has 0(c) = 
4>\ G o a(c + e T ), where e T G G is the identity element, and c G Z>° is arbitrary. 

It remains to combine the two propositions to give an explicit inverse of the isomorphism 
<p\ G (to express the sandpile group G C M C Z>° in terms of configurations), which is part 
(iii) of Theorem 13.41 

Proposition 4.4. The inverse of 4>\ G ^ s given by (13.0.21) . 

Proof. Since <p\ G is an isomorphism and hence a bijection, we need only check that (i) 
the formula (13.0 21) indeed gives a recurrent element (i.e., an element of G) for every v G 

rijeVo ^/ e i.p(i)' and (") ^ ne com P os ition o is the identity. 

The first part (i) is obviously true by definition, using Proposition 14.11 (we defined fj to 
satisfy this). Let us check (ii). To avoid confusion, let us denote the map recursively defined 
in (I3.0.2p by (0 -1 )' for now (so replace every instance of with (0 -1 )'. We then compute 

o (0-7K- = f 3 { Vj ) - J^ir 1 )'^ + ^(0- 1 )'(v), (mod e jiP{j) ) = v r (4.0.11) 

□ 
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